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ANALYTICAL AND EXPERIMENTAL STUDY OF BOILER INSTABILITIES 
DUE TO FEED-SYSTEM - SUBCOOLED REGION COUPLING 
by Thomas M. Grace and Eugene A. Krejsa 
Lewis Research Center 

SUMMARY 

A study was made of instabilities occurring in a forced-flow, electrically heated, 
downflow boiler with constant exit pressure. Steady-state conditions and the frequency 
of oscillation were measured at the onset of oscillations for various feedline character- 
istics. A model representing the dominant dynamic characteristics of the boiler is 
presented. This model was verified by the use of impedance concepts, known character- 
istics of the feed system, and experimental data. 


INTRODUCTION 

Recent emphasis on more efficient and compact methods of heat transfer has led to 
the use of forced convective -boiling systems. One of the more serious problems en- 
countered with these systems is the various types of unstable behavior that occur. 
Instabilities inevitably lead to degraded performance and may become serious enough 
to cause actual failure of the system. It is, therefore, critically important to under- 
stand boiling-loop instability. 

There are a great number of references that discuss experimentally observed boiler 
instabilities. Jeglic and Grace (ref. 1) observed flow oscillations in the subcooled 
region. They found that the onset of these oscillations was strongly dependent on the 
thermodynamic conditions existing in the boiler. In reference 2, Lowdermilk reports 
that oscillations could greatly affect critical heat flux. He found that the instability 
could be prevented by restricting the flow upstream of the test section and that the 
presence of a compressible volume in the feedline had a pronounced destabilizing effect. 
Aladyev and his associates report similar results in reference 3. It has also been 
reported that restricting the boiler exit has a destabilizing effect. Thus, it appears that 
boiler-loop instabilities are caused by an interaction between the feed-system, boiler, 
and exit hardware. 


The dynamic characteristics of the feed-system, boiler, and exit hardware must be 
known in order to describe this interaction. Since little is known about the boiler 
dynamic characteristics, a study was initiated at the NASA Lewis Research Center on a 
boiler facility that has a well-defined feed-system and boiler-exit condition. The facility ' 
was fabricated so that various dynamic characteristics of the feed system could be 
obtained and the pressure at the boiler exit would tend to be constant. Thus, for this 
facility, the instabilities were caused only by feed-system - boiler interaction. 

It is necessary to determine how the feed system and the boiler interact in order to 
make a quantitative analysis of this type of instability. The feed system supplies flow 
to the boiler, and this flow is controlled in part by the pressure at the boiler inlet. Thus, 
the interaction that is a characteristic of this instability is a pressure-flow interaction 
between the feed system and the boiler. Therefore, the characteristics needed to 
analyze the stability of this system are the boiler and feedline impedances. The boiler 
impedance gives the dynamic relation between pressure and flow at the boiler inlet as a 
function of frequency and boiler parameters. The feedline impedance gives the dynamic 
relation between pressure and flow at the boiler inlet as a function of frequency and 
feedline parameters. The feedline impedance can be calculated by using line -dynamics 
theory, but the boiler impedance is unknown. This report presents a model for the 
boiler impedance and verifies the fact that this model represents the dominant charac- 
teristics of the boiler over the range of investigation. 


ANALYSIS OF BOILER DYNAMICS 

If it is assumed that the boiler pressure drop is a function of total flow at the boiler 
inlet and vapor flow at the boiler outlet, then 

< P b>. -< P b> -SP b [w t ,(W T ) ] 

in out out 

(All symbols are defined in appendix A. ) Since (P^) is constant, the impedance of 
the boiler can be written as ou * 
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For a constant -diameter boiler, the vapor generation rate per unit length is 
independent of total flow and depends only on the heat flux and fluid properties. Thus, 
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■ if fluid -property changes are neglected, the vapor flow at the boiler outlet depends only 
on the heat flux and the length of the two-phase region. For electrical heating and thin- 
wall tubing, the heat flux is constant. Thus, changes in the vapor flow at the boiler out- 
let are produced only by changes in the length of the two -phase region: 

(AW ) = 22a ALt 

v out X 

Since the total heated length is constant, changes in the length of the two-phase region 
are produced by changes in the length of the subcooled region: 

iL tp = - aL sc 


Thus, the boiler impedance is given by 
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Solving the subcooled- region energy equation with constant heat flux assumed gives 
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where 


P^SA^sat ~ T j n ) 

ffDq 


is the time required for a particle to travel from the boiler inlet to the end of the sub- 
cooled region, and s is the Laplace variable. (The details of this solution are given 
in appendix B. ) Thus, 
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where 


C„(T 


N — = P 


sat 


T. ) 
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At the neutral stability point s = jco, therefore, the boiler impedance at neutral stability 
is given by 
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This expression describes the change in boiler inlet pressure, caused by a change in 
total flow, as a function of boiler steady-state characteristics and the frequency of 
oscillation. The frequency dependency is contained in the term (1 - e -ja,T )/jwT. The 
amplitude and phase of this term are plotted as a function of cot in figure 1. For the 
data obtained in this study, the term N 36P h /9(W ) | (1 - e”-’ a,T )/jcoT tended to be 
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Figure 1. - Amplitude and phase of frequency dependent term in 
model for boiler inlet impedance. 
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’ motion of the interface between the subcooled and two-phase regions tends to play a 
dominant role. 


APPARATUS AND PROCEDURE 
Description of Apparatus 

The facility used in this study is a closed-loop, low-pressure, forced-flow, 
convective -boiling system. Water flowed vertically down through the boiler. A sche- 
matic diagram of the facility is given in figure 2. The test section (boiler) was con- 
structed of 0. 00167 -foot-thick-wall tubing with a 0. 0208-foot outside diameter. The 
heated length of the test section was 1. 75 feet (100 diam) and was preceded by a 
15-diameter unheated length to reduce the hydrodynamic entrance effects. Heat was 
supplied to the fluid by electrical resistance heating of the tube. Current (ac) was 
supplied to the test section from a saturable reactor controlled by a Variac. A special 
thermocouple was fitted 0. 0208 foot from the exit of the tube as a burnout detector. The 
power to the test section would be shut off automatically if this thermocouple indicated 
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that the wall temperature exceeded a preset value. 

Flow is provided by a gear pump powered by a variable -speed drive. In addition, a 
needle valve was located just downstream of the pump to increase the slope of the 
effective head-flow curve of the pump. This valve was followed by any one of 10 air- 
filled standpipes that acted as an accumulator or hydraulic capacitance. The standpipe 
was teed into the line a short distance downstream of the pump. Following this is a 
0. 04 17 -foot -diameter line leading to a remote controlled globe valve. This valve was 
used to provide variable flow resistance at the inlet to the test section. Provisions 
were made to bleed off the air that would tend to become trapped at the high point in the 
system. 

A condenser-accumulator, a shell and tube heat exchanger, was located immediately 
downstream of the boiler. It served a dual purpose of condensing any vapor produced 
by the boiler and holding the boiler exit pressure constant. The loop was filled in such 
a manner that, when flow was established, a fairly large volume of air was trapped on 
the shell side of the heat exchanger, into which the steam and hot water from the boiler 
were ejected. The water to cool this condenser-accumulator passed through the tubes. 
The air in the condenser held the boiler exit pressure fairly constant. While the boiler 
exit pressure tended to rise as the temperature of the boiler effluent was raised, it rose 
quite slowly and did not change significantly as saturation conditions were approached. 
From all indications, the boiler exit pressure appeared to be insensitive to perturba- 
tions. An additional heat exchanger, with the test fluid circulating through the tubes, 
was used to cool the fluid to a temperature between 60° and 80° F. The temperature 
was determined by the prevailing cooling -tower water temperature. 


Instrumentation 

The instrumentation that was included in this study is shown in the schematic dia- 
gram of the apparatus (fig. 2). Pressure was measured with 0- to 14 400-pound-per- 
square-foot range Bourdon gages. There were four basic steady-state pressure meas- 
urements: Pp just upstream of the globe valve; Pg, at the boiler inlet; Pg, at the 
condenser-accumulator; and P^, at the pump outlet. The pressure drop across com- 
ponents was determined by taking the differences, both in the lines and the system, and 
also accounting for fluid friction in the lines. All the gages had small valves in the 
lines to damp out oscillations and fluctuations, but the valve on the condenser - 
accumulator was always left open. In addition, great care was taken to ensure that the 
lines for Pj, Pg, and P^ were always completely filled with liquid by filling the 
system under vacuum. These lines were then always valved off when flow was not 
established. They would remain shut until flow was circulating in the loop and be 
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closed again before shutting off the pump. This procedure was most critical for gages 
and Pg, which were in the portion of the facility that would fill with air when the rig 
was shut down. A differential gage was used to measure the pressure drop across the 
* pump, but only for ease of operation and as a cross-check on the other measurements. 

A strain -gage -type pressure transducer was mounted in the line about 3 feet downstream 
from the standpipe. It was used to measure pressure oscillations. 

Flow was measured with a turbine flowmeter located just downstream of the pressure 
transducer. The pulsed output of this meter was read directly on a counter to provide 
the basic steady-state flow measurement. In addition, the pulsed output was converted 
to direct current. The converted output was read on a dial meter and recorded on a pen 
recorder. 

Thermocouples were used to measure the temperature of the fluid entering and 
leaving the test section. These temperatures were read on a multipoint recorder. In 
addition, the burnout-detector-thermocouple reading was continuously recorded on a 
stripchart recorder. 

The electrical power supplied to the boiler was measured with a wattmeter. The 
speed at which the pump was rotating was indicated on a dial meter. 

The oscillatory or dynamic data consisted of the flow and pressure measurements 
located slightly downstream of the standpipe. The pressure measurement was taken 
from the strain-gage transducer, and the flow measurement was taken from the output 
of the frequency converter. Both these signals were recorded on a two-channel pen 
recorder. It was possible to amplify both these signals and to eliminate partially the 
direct -current level on the signals. Thus, the oscillatory portion of each signal could be 
played over a wider scale range and, thus, be better examined. A timer fed pulses at 
1 -second intervals to the recorder to provide a time scale for the dynamic data. From 
these measurements, it was possible to determine the onset and the extent of oscillatory 
behavior and the frequency of oscillations. 


Experimental Procedure 

In the course of this program, instability was best approached by increasing the 
power supplied to the boiler. This method was used to obtain the data and gave a fairly 
sharp onset of instability. It was also the method that gave the best control for making 
the small changes needed in approaching the unstable condition. Thus, at any particular 
feedline and mean flow condition, instability was approached by increasing the power 
to the test section. 

Ten different capped-off stainless -steel-tubing standpipes were used as accumula- 

C O 

tors. The standpipes spanned a volume range of 8. 82x10 to 1.62x10 cubic feet. 
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For each accumulator, a series of globe-valve resistance settings were run at three 
different mean flow rates: 1.08x10“^, 2.07x10“^, and 3.06x10"^ slug per second. 

For a given standpipe, the following procedure was used to obtain a data point: 

The pump was turned on to establish flow in the loop, and the valves on gages Pj and 
Pg were opened. Air pockets were removed from the system. The desired mean flow 
rate was then established, and the pump valve was set to take approximately a 5000- 
pounds -per -square -foot drop across it. Then, the globe valve was set to take the 
desired pressure drop across it. At this point, a mean flow rate and a feedline condi- 
tion were established, and power was supplied to the boiler. The power was increased 
in increments and the steady-state data (pressures, power, flow, temperatures, etc.) 
were taken at appropriate intervals. The power was incrementally increased until either 
an oscillation set in or the current limitation on the wattmeter was reached. The onset 
of the oscillations could be observed by examining the pressure and flow signals on the 
pen recorder. The oscillatory signals were recorded for a sufficient length of time to 
establish the frequency of oscillation. Often, however, the oscillation would trigger a 
burnout and automatically shut off the power, sometimes after only a few periods. After 
obtaining a data point in this manner, a new condition would be established and the proce- 
dure repeated. 


DATA REDUCTION 

For a system at neutral stability, the total impedance around the loop must be zero 
at the frequency of oscillation. For this system, the total impedance around the loop is 
the sum of the feedline impedance and the boiler impedance. Thus, at neutral stability, 
the negative of the boiler impedance must equal the feedline impedance. This fact can 
be used to check the validity of the model used to obtain the boiler impedance. 

An expression for the feedline impedance can be obtained by substituting the 
appropriate electrical analogies for the feed-system components. Resistance is 
analogous to pressure -drop devices, such as pumps and valves. Capacitance is 
analogous to compressible volume, and inductance is analogous to fluid inertia. The 
electrical analog of the feed system used in this study is shown in figure 3. The inertia 
of the fluid in the lines from the condenser to the pump is defined as the sum of the 
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Figure 3. - Electrical analog of feed system. 
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ratios of length to area for these lines. The combined resistance of the pump and the 
needle valve just downstream of the pump Rp is the sum of the slopes of pressure drop 
against flow curves for these elements. With air assumed to act as an ideal, isothermal 
* gas, the capacitance of the air in the standpipe C is defined as 


C = 



where V 0 is the initial volume of air at atmospheric pressure P Q , and Pg is the 
mean absolute gas pressure at operating conditions. The inertia of the fluid in the lines 
from the standpipe to the valve at the boiler inlet is defined as the sum of the 
ratios of length to area for these lines. The combined resistance of the valve at the 
boiler inlet and the lines from the standpipe to the boiler inlet R y is obtained from 
slopes of pressure drop against flow curves. The feedline impedance is 


z f = R v + 


^1 - 
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The amplitude and phase of the feedline impedance can be calculated from the following 
relations: 
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fy^p^t^sat “ T iir 

T = 

7rDq 

and where (S6P h /3W t ) is the slope of boiler pressure drop against total flow at 

° (W ) 

out 

constant exit vapor flow, 36P b /3(W ) is the slope of boiler pressure drop 

L out-1 

against exit vapor flow at constant total flow, and u> is the observed frequency of oscil 
lation. The amplitude and phase of the negative of the boiler impedance are obtained 
from the following relations: 

Amplitude = ^ (Real)^ + (Imaginary)^ 



Vapor flow rate at boiler exit, (JN V ) slugs/sec 


Figure 4. - Steady-state boiler pressure drop against vapor flow at boiler exit with constant total flow. 






Figure 5. - Typical plot of impedance in complex plane. Run 794. 
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Phase = - 180° 

V Real / 

Plots of boiler pressure drop against exit vapor flow at constant total flow are 
presented in figure 4. The vapor flow used in these plots was calculated assuming ther- 
mal equilibrium. These plots were used to evaluate the two partial derivatives needed 
to calculate the boiler impedance. All quantities in both the boiler and the feedline 
impedances are evaluated at the observed inceptions of instability. A tabulation of the 
reduced data for each observed inception of instability is given in tables I and II. 

Figure 5 shows typical plots of the boiler and feedline impedances in the complex plane. 

RESULTS AND DISCUSSION 

The feedline impedance, which was determined from known characteristics and 
steady-state data at the inception of instability, is presented in table I for each observed 
point of instability. The negative of the boiler impedance, which was determined from the 
model of boiler dynamic characteristics and steady-state data, is presented in table II 
for each point of instability. Since the sum of the boiler and feedline impedances must 
be zero at neutral stability, the amplitude and phase of the feedline impedance should 
agree with the amplitude and phase of the negative of the boiler impedance if the boiler 
model is valid. 

The amplitude of the boiler impedance is plotted against the amplitude of the feedline 
impedance in figure 6(a). The phase of the negative of the boiler impedance is plotted 
against the phase of the feedline impedance in figure 6(b). The solid lines were drawn at 
45° to facilitate comparison of the two quantities. Most of the amplitudes and phases of 
the boiler impedance agree within ±50 percent of the amplitudes and phases of the feed- 
line impedance. Thus, the data show that the model, while not completely describing the 
boiler, does describe its dominant characteristics. 

While some of the scatter in figure 6 is caused by the evaluation of both the feedline 
and boiler impedances from experimental data, the majority of the scatter is probably 
caused by the assumptions made when the boiler model was derived. Energy storage in 
the tube wall and mass and energy storage in the two-phase region were neglected. 
Although the effect of these assumptions is negligible at low frequencies, it can become 
significant at higher frequencies. In addition to storage, the effects of subcooled boiling 
were neglected. The existence of subcooled boiling should be taken into account when 
the partial derivatives of boiler pressure drop are evaluated. That is, the actual vapor 
flow, rather than the vapor flow calculated assuming thermal equilibrium, should be 
used when evaluating these partials. Subcooled boiling will also affect the change in 
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(b) Phase lag. 

Figure 6. - Concluded. 

vapor flow at the boiler exit due to a change in the two -phase length and will affect the 
value of t. The model used in this report defines r as the time required for a particle 
to travel from the boiler inlet to the point where the bulk fluid temperature is equal to 
saturation temperature. In the presence of subcooled boiling, r should be defined as 
the time required for a particle to travel from the boiler inlet to the point where boiling 
begins. 

CONCLUDING REMARKS 

The analysis of boiler dynamic characteristics presented in this report yields a 
simple model that agrees reasonably well with experimental data. This model shows 
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* that for the conditions of this study, the motion of the interface between the subcooled 
and two -phase regions played a dominant role. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, December 16, 1966, 

120-27-04-27-22. 
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APPENDIX A 


SYMBOLS 


\ 

cross-sectional area of boiler 

P 4 

o 

pump outlet pressure, lb/ft 


tube, ft^ 

q 

heat flux, Btu/(ft 2 )(sec) 

c 

capacitance of air in stand- 

o 

R n 

resistance of pump and needle 


pipe, (ft) (sec ) 

P 

valve downstream of pump, 

C P 

heat capacity of liquid, 


(ft -1 )(sec -1 ) 


Btu/(slug)(°F) 

R v 

resistance of globe valve and 

D 

diameter of boiler tube, ft 


line from standpipe to boiler 

^sc 

length of subcoole4 region in 


inlet, (ft - *)(sec - *) 


boiler, ft 

s 

Laplace variable, sec"*' 

Hp 

length of two-phase region in 

t 

time, sec 


boiler, ft 


liquid temperature, °F 

L 1 

inertia of fluid in line between 
condenser and pump, ft - '*' 

T in 

liquid temperature at boiler 
inlet, °F 

L 2 

inertia of fluid in line between 
standpipe and boiler, ft - * 

T sat 

v ft 

saturation temperature, °F 
initial volume of gas in stand- 

N sc 

C p (T sat - T in>A 


3 

pipe, ft 

(P b). 

m 

pressure at boiler inlet, 
lb/ft 2 

w t 

total flow rate at boiler inlet, 
slug/sec 

( p b ) t 

pressure at boiler outlet, 

(w v ) 

vapor flow rate at boiler outlet, 

° out 

lb/ft 2 

v out 

slugs/sec 

p g 

absolute gas pressure at stand- 

X 

position, ft 

pipe, lb/ft 2 

Z b 

boiler impedance, (ft - *0(sec - *0 

P o 

initial pressure of gas in stand- 
pipe, lb/ft 2 

Z f 

feedline impedance, 
(ft -1 )(sec -1 ) 

P 1 

P 2 

P 3 

pressure just upstream of 
globe valve, lb/ft 2 

pressure at boiler inlet, lb/ft 2 

pressure at condenser- 

6P b 

X 

p l 

boiler pressure drop, lb/ft 2 
heat of vaporization, Btu/slug 

O 

liquid density, slugs /ft 


o 

accumulator, lb/ft 

T 

U) 

subcooled dead time, sec 
frequency, rad/sec 
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APPENDIX B 


DERIVATION OF DYNAMIC CHARACTERISTICS OF SUBCOOLED REGION 

The energy equation for the subcooled region with constant heat flux and constant 
density is 


3T, 


9T 


l 


W t C n — + P Z A t C n — =7rD( l 

1 p ax L * p at 


(Bl) 


or 


7rPq 

ax w t at w t c p 


(B2) 


Putting this equation in perturbation form yields 


8 *r l , 3AT 1 _ -ffDq aw 
0X w t at " CpW 2 1 


(B3) 


Taking the Laplace transform gives 


dAT, P,\s 
+ - - AT 


dX W, 
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7rPq 

C W? 
P t 


AW* 


(B4) 


Integrating gives 


AT^ = Cexp 




„Dq AW t 


C pP l 


s W. 


(B5) 


at X = 0; AT^ = 0. Therefore, 


and 


c . AW t 

w t 


(B6) 
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(B7) 


AT 



-1 


7rDq AW t 


J VA S w t 


Equation (B7) gives the liquid-temperature perturbation at any particular position X as 
a function of the flow perturbation. What is of interest is the change in the subcooled 
length. The change in the subcooled length equals the change in the position at which 
T, = T , . This change, in turn, can be related to the change in temperature at X = L 
through the steady -state temperature profile if second -order effects are neglected: 


AT, 


AL 

sc 


sc 


dT, 


W t C p 

— at 

7rDq 


l 


sc 


dx 


(B8) 


Combining equations (B8) and (B7) gives 



or 


AL 


sc 


sc 


1 - e~ TS \ AW t 

TS / W t 


(BIO) 


where 


p ? A t L sc _ p ^p A t^ T sat ~ T in^ 
W t ?rDq 


(Bll) 


Since 


ffDqLgc ~ ^t^p^sat ~ ^in^ 
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TABLE I. - Continued. AMPLITUDE AND PHASE OF FEEDLINE IMPEDANCE 
AND VALUES NEEDED TO COMPUTE THIS IMPEDANCE 


Run 


Fre- 

quency, 

rad/sec 


Resistance, (ft ^(sec -1 ) 

Pump, 

R P 

Globe value, 
R v 


Inertia of fluid, 

ft" 

1 

■ 

a 


Capacitance 
in stand- 
pipe, 

C, 

(ft) (sec 2 ) 


Feedline impedance 


Amplitude, 

(ft' 1 )(sec' 1 )| 


Phase 

lag, 

deg 


Liquid flow rate at boiler inlet, W t , 2.07XHT 3 slug/sec 


510 

516 

572 

607 

624 


5.71 

5.56 

6.03 

5.15 

4.33 


6.97X10 
7.06 
7.31 
7. 38 
7.47 


2. 32x10 

2.34 

2.28 

2.34 

4.56 


2796 


11 320 


0.83x10 
1.33 
.49 
4.08 
2.28 


-7 


20. 3x10" 
13.3 
30.8 
4.92 
11.1 


66 

69 

61 

58 

58 


650 

656 

715 

727 

734 


2.94 

2.84 

2.79 

3.16 

3.05 


7. 27 
7. 47 
7.23 
7. 10 
7.34 


3.72 

5.80 

5.90 

7.10 

11.00 


♦ ♦ 

3120 33 200 


6.52 
5.31 
3.38 
2.68 

1.52 


6.35 

8.95 

12.0 

13.8 

25.3 


50 

44 

52 

50 

48 


794 

800 

805 

854 

858 


2.62 

2.64 
3.00 

1.65 
1.60 


7. 42 
7.58 
7. 69 
7.27 
7. 34 


5.45 

7.97 

10.76 

3.96 

5.19 


5.93 
4.33 

2.93 
17.14 
14.99 


8.17 
11.9 
16.1 
5. 10 
6.53 


43 

41 

39 

36 

34 


862 

867 

873 

917 

921 


1.65 

1.91 

1.92 
1.26 

.89 


7. 30 
7. 62 
7. 84 
7. 35 
7. 42 


6.85 

9.83 

17.20 

6.99 

9.50 


11.77 

9.30 

4.44 

26.1 

21.75 


8.55 

11.4 

19.4 
7.60 

11.0 


32 

26 

22 

20 

26 


926 

1.68 

7.57 

933 

1.92 

7.72 

973 

00 

CO 

7. 43 

977 

1.14 

7. 39 


13.22 
18. 10 
10.18 
14.00 


♦ 


♦ 


15.3 
8.54 

41.4 
29.0 


13.8 

19.3 

10.9 

14.4 


14 

16 

18 

11 





TABLE H: - AMPLITUDE AND PHASE OF NEGATIVE BOILER IMPEDANCE AND QUANTITIES NEEDED TO 


CALCULATE THIS IMPEDANCE 


Frequency, 

Total heat 

Vapor 

N sc 

86P b 


36P b 


Subcooled 
dead time, 

t , 

sec 

Negative boiler impedance 

rad/sec 

input, 

Btu/sec 

quality 
at boiler 
exit 


aw t 

(tf 1 : 

> 

w . 

out 

Ksec" 1 ) 

9(W V ) 

v out) 
(ft^Hsec 

7 

W t 

- 1 ) 

Amplitude, 

(ft^Hsec- 1 ) 

Phase lag, 
deg 


Liquid flow rate at boiler inlet, W^, 1. 08x10”^ slug/sec 


470 

14.3 

6.46 

0.034 

0.158 

2. 76 X 10 5 

5.8»<10 6 

0.604 

2.71 X 10 5 

137 

471 

13.2 

5.67 


.156 

2.04 

4.32 

.669 

2.22 

141 

473 

4.65 

7.11 

.047 

.159 

3.00 

6.72 

.543 

7.68 

94 

543 

16.0 

7.20 

.054 

. 163 

3.06 

7.20 


3.33 

133 

551 

16.3 

6.95 

.044 

. 162 

3.00 

6.48 

.567 

3.57 

141 

554 

3.21 

5.65 

.004 

.164 

1.68 

3.84 


4.59 

83 

559 

20.0 

9.39 

.104 

.175 

2.58 

10.56 

.453 

4.38 

115 

561 

26.6 

13.54 

.237 

.172 

-8.51 

24.1 

.309 

14.90 

27 

584 

4. 25 

6.58 

.031 

.160 

2.70 

5.64 

.591 

6.52 

117 

598 

4.27 

7.95 

.062 

. 165 

3. 12 

7.80 

.505 

9.46 

79 

690 

2.95 

4.92 

-.006 

. 147 

1.14 

3.12 

.727 

3.38 

79 

695 

2.64 

5.03 

.004 

. 147 

1.68 

3.84 

.710 

4.10 

73 

700 

2.83 

7.47 

mm 

.152 

3.12 

8.28 

.495 

9.46 

52 

762 

2.90 

5.81 

n 

.148 

2.58 

5.16 

.619 

5.41 

73 

767 

2.82 

5.97 

.029 

.148 

2.64 

5.46 

.602 

6.06 

70 

775 

2.86 

8.03 

.084 

.151 


9.24 

.457 

10.80 

47 

779 

3.38 

9.91 

.116 

. 153 

2.16 

11.40 

.408 

14.48 

45 

834 

2.07 

4.98 

-.003 

.151 

1.32 

3.30 

.738 

3.64 

59 

837 

2.06 

4.98 

-.003 

. 152 

1.32 

3.30 

.743 

3.67 

59 

843 

1.90 

6.25 

.037 

.152 

2.82 

6.00 

.591 

6.42 

46 

849 

i 

3.21 

10.17 

.134 

.157 

1.26 

12.6 

.379 

17. 59 

37 

900 

1.40 

4.70 

-.014 

.154 

.60 

2.64 

.798 

3.36 

37 

903 

1.44 

4.72 

-.013 

. 154 

.72 

2.70 

.792 

3.34 

39 

908 

1.51 

6.34 

.028 

. 159 

2.64 

5.46 

.610 

6.12 

37 

957 

1.12 

5.47 

.004 

. 158 

1.68 

3.84 

.702 

4.44 

31 


24 








TABLE II. - Concluded. AMPLITUDE AND PHASE OF NEGATIVE BOILER IMPEDANCE AND QUANTITIES 

NEEDED TO CALCULATE THIS IMPEDANCE 


Run 

Frequency, 

u‘, 

rad/sec 

Total heat 
input, 
Btu/sec 

Vapor 
quality 
at boiler 
exit 

N sc 

36P b 

(Wy) 

out 

(sec" 1 ) 

36P b 

t 

w t 

- 1 ) 

Subcooled 
dead time, 

t, 

sec 

Negative boiler impedance 

3W t 

(ft -1 ) 

3(W y ) 

out 

(ft -1 ) (sec 

Amplitude, 

(ft -I )(sec -1 ) 

Phase lag, 
deg 



Liquid flow rate at boiler inlet, W t , 

3. 06X10"’ 

* slug/sec 



628 

6.16 

13.55 

-0.012 

0.155 

4. 20xl0 5 

4. 32x10® 

0.278 

4. 20X10 5 

94 

634 

6.47 

14. 20 

-.005 

.156 

3.24 

5.63 

.268 

6.16 

73 

640 

5.15 

16. 11 

.008 

.163 

2.58 

5.10 

.246 

4.72 

76 

661 

3.93 

12.20 

-.023 

.152 

6.00 

2.04 

.303 

3.96 

155 

667 

3.88 

13. 20 

-.012 

.152 

4.20 

4.32 

.280 

3.70 

73 

672 

3.55 

15.36 

.006 

.157 

2.70 

7.98 

.249 

9.76 

32 

678 

3.76 

15.09 

.006 

.154 

2.70 

7.98 

.248 

9.50 

34 

686 

4.62 

20.00 

.043 

.157 

6.95 

15.7 

.191 

17.81 

35 

740 

4.14 

14.28 

.005 

. 148 

2.70 

7.85 

.252 

8.88 

39 

745 

4.14 

14.21 

.001 

.148 

2.82 

6.95 

.253 

7.52 

41 

751 

4.33 

15.59 

.018 

.149 

3.12 

10.5 

.232 

12.30 

36 

815 

3.22 

13.80 

.018 

.153 

3.12 

10.5 

.269 

12.78 

31 

819 

3.21 

14. 60 

.002 

.152 

2.86 

7.20 

.253 

8.17 

31 

822 

3.31 

14. 96 

.005 

.153 

2.70 

7.85 

.249 

9.65 

30 

826 

3.98 

18. 90 

.043 

.157 

6.95 

13.1 

.202 

17.88 

32 

829 

3.14 

14.35 

-.005 

.156 

3.24 

5.63 

.264 

5.73 

37 

878 

2.17 

13. 37 

-.013 

.154 

4.20 

4.08 

.280 

2.52 

47 

883 

2.17 

13.80 

-.008 

.153 

3.54 

5.04 

.271 

4.34 

31 

888 

1.89 

14.30 

-.004 

.155 

3.18 

5.94 

.264 

6.07 

22 

892 

3.27 

16.08 

.010 

.158 

2.86 

7.62 

.239 

9.24 

29 

940 

1.50 

13.20 

-.018 

. 157 

4.92 

3.00 

.289 

1.06 

110 

944 

1.65 

13.93 

-.009 

.156 

3.66 

5.76 

.272 

5.40 

22 

948 

2.28 

15.52 

.002 

. 160 

2.86 

7.20 

.251 

8.75 

21 

980 

1.01 

12.52 

-.017 

.150 

4.80 

4.92 

.291 

2.95 

23 
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